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Abstract. Several previous studies highlight pressure (or 
equivalently, pressure altitude) discrepancies between the ra- 
diosonde pressure sensor and that derived from a GPS flown 
with the radiosonde. The offsets vary during the ascent both 
in absolute and percent pressure differences. To investigate 
this problem further, a total of 731 radiosonde/ozonesonde 
launches from the Southern Hemisphere subtropics to north- 
ern mid-latitudes are considered, with launches between 
2005 and 2013 from both longer term and campaign-based 
intensive stations. Five series of radiosondes from two man- 
ufacturers (International Met Systems: iMet, iMet-P, iMet-S, 
and Vaisala: RS80-15N and RS92-SGP) are analyzed to de- 
termine the magnitude of the pressure offset. Additionally, 
electrochemical concentration cell (ECC) ozonesondes from 
three manufacturers (Science Pump Corporation; SPC and 
ENSCI/Droplet Measurement Technologies; DMT) are ana- 
lyzed to quantify the effects these offsets have on the cal- 
culation of ECC ozone (O 3 ) mixing ratio profiles (O 3 MR) 
from the ozonesonde-measured partial pressure. Approxi- 
mately half of all offsets are > ±0.6 hPa in the free tropo- 
sphere, with nearly a third > ±1.0 hPa at 26 km, where the 
1.0 hPa error represents ~ 5 % of the total atmospheric pres- 
sure. Pressure offsets have negligible effects on O 3 MR be- 
low 20 km (96 % of launches lie within ±5 % O 3 MR error at 
20km). Ozone mixing ratio errors above lOhPa (~30km), 
can approach greater than ±10% (>25% of launches that 
reach 30 km exceed this threshold). These errors cause dis- 
agreement between the integrated ozonesonde-only column 


O 3 from the GPS and radiosonde pressure profile by an av- 
erage of ±6.5 DU. Comparisons of total column O 3 between 
the GPS and radiosonde pressure profiles yield average dif- 
ferences of ±1.1 DU when the O 3 is integrated to burst with 
addition of the McPeters and Labow (2012) above-burst O 3 
column climatology. Total column differences are reduced to 
an average of —0.5 DU when the O 3 profile is integrated to 
10 hPa with subsequent addition of the O 3 climatology above 
10 hPa. The RS92 radiosondes are superior in performance 
compared to other radiosondes, with average 26 km errors 
of — 0.12hPa or ±0.61% O 3 MR error. iMet-P radiosondes 
had average 26 km errors of — 1.95hPa or ±8.75% O 3 MR 
error. Based on our analysis, we suggest that ozoneson- 
des always be coupled with a GPS-enabled radiosonde and 
that pressure -dependent variables, such as O 3 MR, be recal- 
culated/reprocessed using the GPS-measured altitude, espe- 
cially when 26 km pressure offsets exceed ±1.0hPa/±5 %. 


1 Introduction 

A number of fundamental intercomparison studies of ra- 

diosonde (e.g., Nash et ah, 2006, 2011; da Silveira et ah, 

2006) and ozonesonde (e.g., Smit et ah, 2007; Deshler et 

ah, 2008) instrument performance have appeared within the 
past two decades. Radiosonde investigations have focused 

on comparisons of instrument type with respect to tem- 
perature (Gaffen, 1994; Gaffen et ah, 1999; Steinbrecht et 
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al., 2008; Sun et al., 2010), humidity (Vomel et al., 2007; 
Yoneyama et al., 2008; Miloshevich et al., 2006; Sun et al., 
2010) and pressure (De Muer and De Backer, 1992; Inai 
et al., 2009; Hurst et al., 2011) measurements and typically 
have been associated with the adoption of new sonde mod- 
els. The performance of electrochemical concentration cell 
(ECC) ozonesonde instruments, of which there have been 
three manufacturers since the 1970s, has been compared with 
various compositions of sensing solution type in laboratory 
conditions (Smit and Kley, 1998; Smit et al., 2007; Smit and 
Berg, 2011), and field conditions (Komhyr et al., 1995a, b; 
Thompson et al., 2007; Deshler et al., 2008). The discrep- 
ancies among the ozonesonde instrument-sensing-solution 
combinations are ~ 5 — 15 % relative to an absolute O 3 mea- 
surement, depending on ECC manufacturer, and are pressure 
(and thus, altitude) dependent. The O 3 community has made 
many attempts to homogenize standard operational proce- 
dures (Deshler, 2012; WMO, 2013) for station pre-flight 
preparations and intercomparison of different ECC cells, so 
some of the performance characteristics of the ozonesonde 
prior to launch are well understood. At present, the global 
ozonesonde community is reprocessing thousands of O 3 pro- 
files from dozens of stations to produce a more accurate pro- 
file data set for trend analysis (Ozonesonde Data Quality As- 
sessment, 03S-DQA; Smit et al., 2012). In this effort, the 
pressure measured by the radiosonde to which the O 3 partial 
pressure is referenced, has been taken as free of biases. 

The relatively recent widespread use of GPS-enabled ra- 
diosondes has shown that pressure sensors often differ from 
the pressure derived from the GPS data. These errors propa- 
gate to errors in the calculated O 3 mixing ratio (C>3mr)- 

1.1 Efforts to quantify radiosonde errors and biases 

Numerous intercomparison studies investigate biases in the 
pressure, temperature, humidity and GPS measurements 
amongst various radiosonde types, da Silveira et al. (2006) 
launched five types of GPS-enabled radiosondes in groups 
to analyze GPS measurements in addition to meteorological 
measurements. They found the reproducibility and compar- 
isons of GPS altitude in the stratosphere were within ±20 m. 
Similar results were obtained by Nash et al. (2006, 2011), 
who, in two different studies, found pressure sensors to be su- 
perfluous based on excellent performance of GPS radioson- 
des. Hurst et al. (2011) compared RS92 and iMet pressure 
measurements and found that paired RS92 radiosondes all 
compared to within ±0.3 hPa in the stratosphere and that 
iMet radiosondes averaged approximately 0.8 hPa lower than 
the RS92s between 25-30 km, an error of >5%. Inai et 
al. (2009) studied individual RS80 radiosonde launches to 
compare pressure derived from GPS measurements with the 
radiosonde pressure sensor and found pressure sensor biases 
of — 0.5hPa above 20 km. These pressure errors need to be 
considered in the context of O 3 MR measurements and total 
column O 3 integration. 


Lately, radiosonde manufacturers (e.g., Lockheed Martin 
Sippican, Inc., GPS Mark II Microsonde) have been produc- 
ing radiosondes without pressure sensors, relying on GPS al- 
titude, temperature, and humidity measurements and the hyp- 
sometric equation to derive pressure data. This same tech- 
nique is used in this study and will be described below. 

1.2 Importance of accurate O 3 measurements 

The importance of long-tenn, accurate O 3 profile records is 
well documented in climate reports (IPCC, 2007), O 3 assess- 
ment reports (WMO, 2011), and numerous studies of trends 
in tropospheric (Logan et al., 1994; Logan et al., 1999, 2012; 
IPCC, 2007), stratospheric (Miller et al., 1995; Froidevaux et 
al., 1996; Liu et al., 2006; Rault and Taha, 2007; Jiang et al., 
2007; Kroon et al., 2011) and total column O 3 (Thompson 
et al., 2003, Osterman et al., 2008). Furthermore, ozoneson- 
des provide the highest vertical resolution (~ 100 m) O 3 
measurements from the surface to over 30 km. For this rea- 
son, the satellite remote-sensing community continues to use 
ozonesonde profile data for validation and improvement of 
O 3 profile retrievals (e.g., Nalli et al., 2013). The absolute ac- 
curacy of radiosonde measured pressure profiles themselves 
also has ramifications in the validation of satellite-derived 
pressure-profile Environmental Data Records (EDRs; Nalli 
et al., 2013). 

Biases in O 3 measurements from the use of several differ- 
ent types of ECC ozonesonde manufacturers, as well as dif- 
ferent potassium iodide sensing solution strengths and sonde 
preparation techniques have made the homogenization of the 
historical ozonesonde record a necessity. The goals of the 
homogenization process performed through collaborative ac- 
tivities from WMO’s 03S-DQA (Smit et al., 2012) and SI2N 
(this special issue) are to compile the highest accuracy O 3 
profile records for more robust trend studies and satellite 
comparisons (Deshler, 2012). With the ongoing reprocess- 
ing of ozonesonde data, it is vital to identify every potential 
bias or error in the O 3 measurements. A goal of this paper 
is to contribute to a consensus-based recommendation on the 
handling of these radiosonde errors. Note, however that the 
data we present here and recommendations made apply only 
to radiosondes launched in the GPS era. 

In the present investigation a series of 624 ozonesonde - 
radiosonde instrument packages and 107 RS92-SGP ra- 
diosondes flown solo have been analyzed. In this paper, we 
address the following questions: 

1. What are the statistical characteristics for pressure dif- 
ferences (“offsets”) between the pressure sensor and 
that derived from the GPS? How do the offsets vary as 
a function of pressure (altitude)? 

2. How do the offsets vary between radiosonde models? 
In this study, we analyze the RS92-SGP, three versions 
of International Met Systems (iMet) radiosondes, all of 
which have GPS units integrated into the radiosonde 
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hardware, and the RS80-15N flown with a separate 
Garmin GPS unit attached inside the ozonesonde sty- 
rofoam box. 

3. In addition to pressure offsets, some of the radiosondes 
demonstrate highly variable pressure measurements 
during ascent, especially in the stratosphere. What are 
the statistical characteristics of this variability? 


(Steinbrecht et al., 2008), to a more accurate solid-state sil- 
icon barocap sensor. Note that the quoted Vaisala pressure 
accuracies are valid only if a ground check with an indepen- 
dent surface pressure measurement and calibration are per- 
formed prior to launch. The RS92 radiosondes used in this 
study underwent this check. 

2.2 Ozonesonde measurements 


4. How do the radiosonde pressure offsets propagate to 
the O 3 profiles? How is integrated total O 3 to either 
the balloon burst altitude or a pressure cut off (e.g., 
~ 11 hPa/~30km, as recommended in Dobson, 1973 
or lOhPa, as utilized in Thompson et al., 2003, 2007), 
and an extrapolated add-on determined from a clima- 
tology like McPeters and Labow (2012) affected? 

The soundings were taken in the 2005-2013 period in a 
range of locations from the northern mid-latitudes through 
the subtropics and tropics to southern subtropics (Table 1). 


2 Methodology 

2.1 Site and instrument descriptions 

A total of 73 1 radiosondes were analyzed for this study, with 
ozonesonde/radiosonde pairs accounting for 624 of those 
profiles. Our analysis includes data from twelve different 
launch sites (including two simultaneously operated, closely 
located sites in Houston, TX) launching five types of ra- 
diosondes, and spanning the years 2005-2013 (Table 1). The 
locations range from the southern subtropics (Irene, South 
Africa) to the northern mid-latitudes (Sapporo, Japan) with 
every month of the year represented. Stations include both 
those making regular ozonesonde launches (Irene, Houston, 
Beltsville) and those making intensive launches for specific 
campaigns (see Table 2 for campaign details), as well as other 
profiling missions at other sites. 

Two radiosonde types from Vaisala (Vantaa, Finland; 
RS80-15N, RS92-SGP; herein RS80 and RS92) and three 
from International Met Systems (Grand Rapids, MI, USA; 
iMet, iMet-P, iMet-S) were launched at the various loca- 
tions. Analyses are presented for each radiosonde type. 
The number of launches of each radiosonde type and the 
manufacturer-quoted pressure accuracies/uncertainties are 
given in Table 3 . International Met Systems uses a piezoresis- 
tive silicon device to measure pressure and quotes only one 
pressure accuracy throughout the manufacturing of their ra- 
diosondes from 2009-2013. The analyses are still presented 
by each series type (based on serial numbers that are, in 
general, temporally partitioned) to determine any differences 
throughout the evolution of iMet radiosonde production. The 
RS92 radiosondes received a significant pressure sensor up- 
grade from the RS80s, moving from an aneroid capacitor, 
which is observed to have a low bias in the stratosphere 


Each of the Science Pump Corporation (SPC) and EN- 
SCl/Droplet Measurement Technologies (DMT) ozoneson- 
des in this study operate using the electrochemical concen- 
tration cell (ECC; Komhyr, 1969) technique where ambient 
air is bubbled through a potassium iodide solution. The sub- 
sequent reactions generate two electrons per O 3 molecule, 
so the current measured through an attached circuit board is 
converted to O 3 partial pressure (/ 703 ) via this equation: 

P 03 = 0.043085 ■ (I M - /b), (1) 

(e x F) 

where the constant 0.043085 is derived from the ratio be- 
tween the gas constant R and the Faraday constant, 7 p the 
measured temperature of the ozonesonde pump, e the pump 
efficiency, F the volumetric flow rate through the pump, 7 m 
the measured electrical cell current, and 7 b the background 
cell current quantified in lab testing prior to launch (WMO, 
2013). 

Since O 3 MR is calculated from po 3 and total air pressure, 
Pair • 

03MR = 7 ^ 03 / Pair, (2) 

any bias or error in the radiosonde pressure measurement in- 
troduces error in C>3mr- The /?o3 measurements have typ- 
ical tropospheric accuracies on the order of —7 to +17%, 
improving to ±5 % in the low to mid-stratosphere with de- 
creasing accuracy above lOhPa, provided standardized and 
accepted ozonesonde conditioning and launch procedures are 
followed (Komhyr et al., 1995b; WMO, 2013). 

2.3 Calculation of GPS pressure 


The pressure altitude reported by the radiosonde is given 
in geopotential altitude (Z), using standard gravity (go = 
9.80665 ms -2 ). Conversely, the GPS altitude is reported 
as a geometric altitude (77), and the latitude (g^ » 9.78- 
9.83 ms -2 ) and altitude-dependent gravity is used to calcu- 
late pressure. The equation for gravity with latitude (</>) is 
estimated from the WGS-84 ellipsoid (National Imagery and 
Mapping Agency, 2000, p. 42): 


1 + 0.00193185138639sin 2 (/> 
g 0 = 9.7803267714 • == 

y 1 - 0.00669437999013 sin 2 </> 


The gravity with altitude is then given by 


gc/>,H = g<j> + [ 


Gm^ 

(r E + H ) 2 


Gm^ 



], 


(3) 

(4) 
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Table 1. Balloon launch locations with latitude/longitude coordinates, number of launches, radiosonde types used and lengths of records 
used in this study. 


Location 

Lat/Lon 

Launches 

Radiosonde types 

Length of record 

Irene, South Africa 

— 25.91°/28.21° 

28 

RS92-SGP 

20 Feb 2012-8 May 2013 

Las Tablas, Panama 

7.75°/— 80.25° 

23 

RS80-15N 

13 Jul 2007-9 Aug 2007 

Houston, Texas (Two Locations) 

29.72°/— 95.34° and 
30.03°/— 94.08° 

275 

RS80-15N, iMet, iMet-P, iMet-S 

20 Sep 2005-26 Jan 2013 

Ronald H. Brown R/V, Gulf of 
Mexico 

24.8° to 29.7°/— 94.7° 
to -83.5° 

37, 107 

radiosonde 

only 

RS92-SGP 

27 Jul 2006-1 1 Sep 2006 

Ronald H. Brown R/V, AEROSE 
Expeditions 

—23.5° to 31.8°/— 76.0° 
to -5.9° 

69 

RS92-SGP 

1 1 May 2007-22 May 2010 

Idabel, Oklahoma 

33.89°/— 94.75° 

57 

iMet, iMet-P, iMet-S 

19 Jul 2010-6 Oct 2012 

Porterville, California 

36.03°/- 119.05° 

25 

iMet, iMet-S 

16 Jan 2013-6 Feb 2013 

Beltsville, Maryland 

39.05°/— 76.88° 

16 

RS92-SGP 

27 Jun 2007-7 Aug 2007 

Edgewood, Maryland 

39.41°/— 76.30° 

36 

iMet-S 

28 Jun 201 1-30 Jul 2011 

Valparaiso, Indiana 

41.46°/— 87.04° 

31 

RS80-15N 

19 Apr 2006-3 Nov 2007 

Sapporo, Japan 

43.07°/141.35° 

27 

RS80-15N 

6 Aug 2008^4 Sep 2009 


Table 2. List of campaigns, their respective locations, websites and radiosonde types launched. 


Campaign 

Location(s) 

Website 

Radiosonde 

types 

Dates available 

TC 4 : Tropical Composition, Cloud, and 
Climate Coupling 

Las Tablas, Panama 

http://www.nasa.gov/ 

mission_pages/TC4/ 

RS80-15N 

Jul-Aug 2007 

IONS-06: INTEX-B Ozonesonde Network 
Study 2006 

Houston, TX, Ronald H. Brown 
R/V 

http://croc.gsfc.nasa. 

gov/intexb/ions06.html 

RS80-15N, 

RS92-SGP 

Mar-May, Aug-Sep 2006 

DISCOVER- AQ: Deriving Information on Sur- 
face Conditions from Column and Vertically 
Resolved Observations Relevant to Air Quality 

Edgewood, MD, Porterville, 
CA 

http://www.nasa. 
gov/mission_pages/ 
discover- aq/index . html 

iMet, iMet-S 

Jul 2011, Jan-Feb 2013 

AEROSE: NOAA Aerosols and Ocean Science 
Expeditions 

Ronald H. Brown R/V 

http://aerose.org/ 

RS92-SGP 

May 2007-May 2010 

TOPP: Tropospheric Ozone Pollution Project 

Houston, TX, Idabel, OK, Val- 
paraiso, IN, Sapporo, Japan 

http://physics.valpo. 

edu/ozone/index.html 

RS80-15N, 
iMet, iMet-P, 
iMet-S 

Sep 2005-Jan 2013 


where g,/, is the surface gravity at a given location, G the 
gravitational constant (6.67428 x 10“ 11 N m 2 kg -2 ), the 
mass of Earth (5.97 3 6 x 10 24 kg), ry the average radius of 
Earth (6.371 x 10 6 m), and H the GPS altitude from the ra- 
diosonde. 

This process is the reverse of obtaining a geopotential al- 
titude from the radiosonde pressure measurements, but with 
a geometric altitude. We note that the reported GPS altitude 
is actually an ellipsoidal altitude, though the difference be- 
tween that and altitude AMSL (geoidal altitude; National Im- 
agery and Mapping Agency, 2000, p. 68) is reconciled with 
the input of the station AMSL altitude as the initial GPS al- 
titude prior to launch. Surface pressure from the radiosonde 
(often set at the launch site from a high-precision barometer) 
is used to initialize the GPS pressure calculation from the 
hypsometric equation: 


Pops/ = PGPS/-iexp[— 


g<t>.n AH 

Rd Tv avg 


( 5 ) 


Here, pops* is the pressure calculated from g,/,.//. the lat- 
itude and altitude-dependent gravity, A H, the change in ge- 
ometric GPS altitude from consecutive measurements, P c |, 


the specific gas constant for dry air (287.05 J kg -1 K -1 ), 
and r vav g, the average virtual temperature of the consec- 
utive measurements. Calculating pressure in iterative fash- 
ion from measurement to measurement throughout the pro- 
file reduces the error that use of a standard atmosphere or 
scale height would introduce. Since the uncertainty in the 
GPS altitude is small (Nash et al., 2006, 2011), usually 
within ±20 m (Vaisala RS92 technical specifications testing) 
to ±30 m (iMet radiosonde 2a error specifications), the un- 
certainty resulting from GPS altitude measurements in the 
calculated /?gps will be negligible in the stratosphere. An 
additional source of uncertainty in pops results from errors 
and biases in radiosonde temperature and humidity measure- 
ments (Richner and Viatte, 1995; Hurst et al., 2011). Large 
systematic biases in radiosonde temperature measurements 
can cause some errors in the calculated pgps profile. How- 
ever, the characteristic pressure errors resulting from temper- 
ature biases bear no resemblance to the errors seen in this 
paper - therefore we rule out this factor as the cause of the 
pressure offset. The pops calculation assumes that the atmo- 
sphere is in hydrostatic balance with pressure dependence 
only in the vertical and with negligible changes horizontally 
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Fig. 1. Edgewood, MD iMet-S profile from 14 July, 2011 of GPS 
and pressure altitude differences (black), and pressure differences 
after recalculation of pressure data from GPS measurements (grey). 
The red dashed line marks the zero line for reference. The median 
26 km pressure offset for this launch is —1.17 hPa. 


Fig. 2. Sapporo, Japan, RS80 profile from 19 August 2008 showing 
original pressure (blue) and recalculated GPS pressure (red) O 3 MR 
profiles. The inset figure is the same profile, zoomed in to highlight 
O 3 differences in the stratosphere. The median 26 km pressure off- 
set for this launch is — 1 .94 hPa. 


and in time. The radiosonde makes the same assumption 
when deriving a geopotenlial altitude from the pressure mea- 
surements through use of the hypsometric equation. 

An example of the differences in radiosonde pressure 
(herein p) and pops (treated as the reference), as well as 
the pressure altitude and GPS altitude differences, are shown 
in Fig. 1. Large differences, on the order of hundreds (and 
sometimes thousands) of meters, between pressure altitude 
and GPS altitude are an indication of systematic errors in re- 
ported pressures. For the remainder of this paper, we define 
the pressure offset to be p-pGvs- 

Variability in the pressure offset appears in the lower tro- 
posphere since a difference of just a few meters between the 
GPS and the pressure altitude can cause several tenths of 
1 hPa difference between the calculated and measured pres- 
sures. The noise in the pressure offset stabilizes in the strato- 
sphere and tends to remain somewhat constant until balloon 
burst. 

2.4 Recalculation of pressure-dependent data 

Using the pressure calculated from the GPS measurements, 
any pressure-dependent variables can be recalculated and 
compared to the original measurements. In addition to the 
reported altitude and pressure differences between the GPS 
and radiosonde measurements, we examine the effects on the 
O3MR and total column O3. (Note that the pressure correc- 
tions implemented here also result in a need to recalculate po- 
tential temperature and, to a lesser extent, water vapor mix- 
ing ratio, but we do not discuss these modifications here). We 
choose to examine O3MR rather than /?o3 so we can describe 
statistics on the differences between “coincident” measure- 
ments on a single GPS altitude, that is, the original O3MR and 
the GPS-pressure derived and recalculated C>3 mr- The po3 is 


not dependent on ambient pressure (except for minor pump 
efficiency corrections) so coincident measurements will not 
change in magnitude, and only an altitude shift in the profile 
will be evident. 

The recalculation of O3MR causes differences in both the 
O3 magnitude and profile shape, particularly above 26 km 
and near the burst altitude (Fig. 2). Depending on the severity 
of the pressure offset, O3MR errors can approach ±1-2 ppmv 
(parts per million by volume; ±10-20 % error) or greater in 
the stratosphere in the heart of the ozone layer. Differences 
between GPS altitude and pressure altitude can cause the ap- 
parent O3MR maximum to shift by as much as ±2 km, hav- 
ing further consequences for stratospheric satellite measure- 
ments and comparison/validation studies with ozonesondes. 

We note that a pressure-dependent pump correction fac- 
tor (PCF) is applied to /?o3 based on decreasing ozonesonde 
pump efficiency in the stratosphere, particularly above 25 km 
(Johnson et al., 2002). Flowever, both the application of var- 
ious PCFs in different processing software and the negligi- 
ble (~ 0.5 % difference in PCF between 20 and 18 hPa, near 
where statistics from this paper are presented) difference the 
PCF has between p and pops profiles lead us to neglect this 
small correction. This effort to quantify pressure sensor er- 
rors is separate from and not a substitute for the PCF prob- 
lem. 

3 Results 

The IONS-06 campaign (Table 2) in March-May and 
August-September 2006 provided an opportunity to compare 
coincident O3 profiles from the University of Houston Main 
Campus (UH) and the Ronald H. Brown (herein RHB), op- 
erated by NOAA to record profiles near the Houston Ship 
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Fig. 3. Nearly coincident profiles from 30 August 2006 from the 
Ronald H. Brown (RS92), and Houston, TX (RS80). The original 
C>3MR profiles are shown in black (Houston) and red (RHB). The 
inset highlights improved stratospheric O 3 MR agreement from the 
coincident RHB and Houston profiles after GPS reprocessing with 
the corrected profiles from pqps shown in grey (Houston) and or- 
ange (RHB). The median 26 km pressure offsets for these launches 
are —1.55 hPa (Houston) and 0.00 hPa (RHB). 


Channel and in Galveston Bay. The comparisons allow us to 
test confidence in the paps recalculation procedure, namely 
the reproducibility of stratospheric O 3 MR using ozoneson- 
des with different radiosonde types released closely in space 
and time. Nine such pairs occurred within 90 min of each 
other in IONS-06, with RHB launching RS92s and the UH 
site launching RS80s with a separate GPS unit attached. An 
example of one pair, 15 min and 77 km apart on 30 Au- 
gust 2006, is shown in Fig. 3. The two profiles show simi- 
lar tropospheric O 3 MR with or without correcting the pres- 
sure offset (the p and pops profiles are indistinguishable be- 
low ~ 15 km). The GPS corrected pressure, however, results 
in better agreement in stratospheric C> 3 mr- Before correcting 
the pressure offset error, the mixing ratio differences between 
the two flights are greater than 1 ppmv near the UH balloon 
burst altitude (also note the altitude shift; Fig. 3). Those dif- 
ferences become markedly smaller to within 0. 1-0.2 ppmv 
after correction of both profiles using pops- Both the shift 
in the altitude and correction of the O 3 MR contribute to this 
improved agreement. 

3.1 Statistical characteristics of the pressure offsets 

The median pressure offset for each km altitude bin (as in 
Hurst et al., 2011) from 1-30 km is shown in Fig. 4. The 
tight grouping of RS92 launches about the zero line is distin- 
guishable, with considerably more spread near the top of the 
profiles measured with the other radiosonde types. Most indi- 
vidual radiosondes show less variable pressure offsets in the 


stratosphere, with the RS92s converging to zero. The iMet-P 
radiosondes exhibit a peculiar S-shape pressure offset peak 
around 5 km that is not understood (we can find no artifact or 
geophysical cause). 

At 26 km (an altitude 69% of profiles reach, also chosen 
because p ~ 20 hPa at 26 km), the iMet and RS80 radioson- 
des exhibit the most variable pressure offsets, with mean 
offsets of —0.93 hPa and — l.OlhPa, respectively (see Ta- 
ble 4 for offset averages and percentiles). In Fig. 5, we see 
the radiosonde-measured pressure is consistently lower than 
POPS for many of the radiosonde types, and nearly a third of 
all launches have an offset of ±1.0hPa at 26 km. The least 
variability is exhibited by the RS92s with only a —0.12 hPa 
average offset and just two outlier profiles beyond ±1.0 hPa 
at 26 km. 

Figure 6 shows pressure offsets at various altitudes as a 
function of the pressure offset at the burst altitude. The vari- 
ance within the figure at different altitudes implies that the 
pressure offsets are not constant throughout most of the pro- 
file, and that a constant pressure correction cannot be ap- 
plied to the entire profile. Only when the balloon reaches 
the stratosphere and around 15-20 km is a strong relationship 
evident. The tropospheric offsets appear much less constant 
than the stratospheric offsets, likely from variability in the 
GPS altitude and pressure sensor causing significant noise in 
the pressure offset below 10 km. In the troposphere an offset 
of several hPa only represents a few percent of the total at- 
mospheric pressure. At 20 km, 96 % of all launches have less 
than ±5 % error. As a result, the true magnitude of the pres- 
sure offset and need for reprocessing cannot be determined 
until the balloon has reached the stratosphere (see Fig. Al in 
Appendix A for altitude differences with pressure offset). 

3.2 O 3 MR offsets 

Pressure offsets of only a few tenths of 1 hPa are the equiv- 
alent of 5-10% errors in the total atmospheric pressure at 
the balloon burst altitude near 30 km. This pressure offset er- 
ror results in an error in the calculated O 3 MR of the same 
magnitude (Fig. 7). We define the O 3 MR offsets as [C> 3 mr(/j)- 
03mr(GPS)]/C>3mr(GPS)- Figure 7 demonstrates how a nearly 
constant stratospheric pressure offset results in an O 3 MR off- 
set that grows in magnitude with altitude, with many profiles 
beyond ±10% error in the stratosphere. At such magnitudes, 
this error becomes a significant component of the overall er- 
ror budget associated with O 3 profile data from ozonesondes, 
and is beyond the intrinsic uncertainty of the O 3 measure- 
ments. 

Table 4 examines the O 3 MR errors by radiosonde type. As 
with the pressure offsets, the most variable O 3 MR percent off- 
sets are displayed by the iMet and RS80 radiosondes with 
±4.42% and ±4.75%, at 26 km, respectively. The iMet-P 
launches have an average offset at 26 km of ±8.75% that 
increases to ±15.9% by 30 km, leading to an average error 
greater than 1 ppmv O 3 MR by balloon burst. This large error 
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Fig. 4. Median pressure offset (p-pcj>s) for every 1 km altitude bin from 1-30 km for each radiosonde type (grey). Average offsets (black 
solid line) for each grouping of radiosondes are shown along with 10th and 90th percentiles (black dashes). A red dashed line marks the zero 
line for reference. 


Table 3. Radiosonde types with number of launches, quoted pressure uncertainties/accuracies from the manufacturer, and dates of available 
launches. We note the various iMet series have had no appreciable changes to the pressure sensors, but are split in these analyses for 
convenience and ease of interpretation. 


Radiosonde type 

Launches 

Quoted pressure uncertainty/accuracy 

Length of record 

iMet 

106 

1070-400 hPa: 1.8 hPa/400-4 hPa: 0.5 hPa 1 

28 May 2009-6 Feb 2013 

iMet-P 

52 

1070-400 hPa: 1.8 hPa/400^1 hPa: 0.5 hPa 1 

23 Mar 2012-26 Jan 2013 

iMet-S 

69 

1070-400 hPa: 1.8 hPa/400-4 hPa: 0.5 hPa 1 

5 Nov 2009-16 Jan 2013 

RS80-15N 

247 

1080-3 hPa: 1.0 hPa 2 

20 Sep 2005-23 Apr 2011 

RS92-SGP 

257 

1080-100 hPa: 0.5 hPa/100-3 hPa: 0.3 hPa 2 

27 Jul 2006-8 May 2013 


* The iMet values given are 2 a accuracy limits. - The RS80 and RS92 values given are 2 a limits on sounding reproducibility, valid only after 
performing a ground check between the radiosonde and an independent measurement of surface pressure. 


at 30 km is common; over a quarter of all launches that reach 
this altitude have O3MR errors > ±10%. 

Two distinct offset regimes are detected in the RS92s in 
Figs. 4 and 7, separable mainly by the launch sites Beltsville 
(one summer of data) and RHB near Galveston Bay (single 
campaign; see Figs. A2 and A3 in Appendix A for pres- 
sure and O3MR offsets by launch site). The Beltsville pres- 
sure offsets lie slightly to the left of the zero line, and the 
RHB offsets straddle the zero line. Similar offset group- 
ings are also observed in the campaign-based launches from 
Porterville, CA (iMet, only one iMet-S), Las Tablas, Panama 
(RS80) and the set of iMet-P sondes launched in the course 
of 10 months at Idabel and Houston. This suggests that 
particular “batches” of radiosondes, regardless of manufac- 
turer/type, may have offsets that generally behave in similar 


manners. As a result, we caution against drawing conclusions 
about radiosonde types (particularly iMet-P radiosondes in 
this study) from offsets appearing in only one set or batch of 
sondes. 

3.3 Column ozone measurements 

Because the pressure offset affects both the apparent al- 
titude and magnitude of O3MR, it is also of interest to 
compute the influence on total column amount of O3. 
Each ozonesonde that reached 26 km was integrated to ob- 
tain a column O3 amount in Dobson Units (1 DU = 2.69 x 
10 16 molecules cm -2 ) from both the original pressure profile 
and the recalculated paps profile. As expected, considerable 
differences in the column integrated to the sonde burst alti- 
tude appear closely related to the pressure offset magnitude 
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Table 4. Various pressure and O 3 statistics separated by radiosonde type. All columns are presented in 10th percentile, mean and 90th 
percentile values. Values are reported as original pressure-profile data minus GPS-calculated pressure-profile data. 


Radiosonde 

type 

Pressure 
offset 
(hPa, 26 km) 

O 3 MR 
error 
(%, 26 km) 

Sonde column 
difference 
(to burst, DU) 

Sonde column 
difference 
(to lOhPa, DU) 

Total column 
difference 
(to burst + add-on, DU) 

Total column 
difference 
(to 10 hPa + add-on, DU) 

iMet 

-1.85, -0.93,0.22 

-0.88,4.42,8.94 

-0.1, 11.2, 19.4 

- 0 . 1 , 3 . 5 , 7.0 

- 1 . 0 , 3 . 8 , 7.8 

- 1 . 4 , 0 . 8 , 3.5 

iMet-P 

-2.86,-1.95,-1.00 

4.47,8.75, 13.3 

5.5, 11.8, 16.8 

5.5, 11.6, 16.3 

- 3 . 6 , - 1 . 6 , 0.7 

-3.6, -1.7, 0.7 

iMet-S 

-1.36, -0.58,0.46 

-1.80,2.49,5.65 

-3.2, 4.9, 17.0 

-3.0, -0.1, 10.6 

- 1 . 1 , 1 . 7 , 3.2 

-0.9, -0.1, 2.6 

RS80-15N 

-2.49, -1.01,0.03 

-0.13,4.75, 11.0 

-0.2,9.3,24.6 

-1.2, 3.5, 11.7 

- 1 . 5 , 1 . 9 , 6.1 

-1.9, -0.6, 1.8 

RS92-SGP 

-0.47, -0.12,0.07 

-0.25,0.61, 1.99 

-1.5, 0.5, 2.3 

-1.6, -0.6, 0.4 

- 1 . 1 , - 0 . 2 , 0.8 

-1.6, -0.6, 0.3 

All Sondes 

-1.99, -0.72,0.08 

-0.33, 3.56,9.56 

-1.1, 6.5, 16.6 

-1.5, 2.6, 11.5 

- 1 . 6 , 1 . 1 , 3.9 

-1.9, -0.5, 1.8 



Fig. 5. Histogram of 26 km pressure offset in percent frequency 
hy radiosonde type. Data are binned every 0.5 hPa. The various ra- 
diosonde types are identified by their respective colors. 


in the stratosphere (Fig. 8a) - the radiosonde types that dis- 
played the largest pressure and O 3 MR offsets also present the 
largest sonde column offsets. The iMet, iMet-P and RS80 
sonde-only column O 3 is consistently ~ 10 DU higher than 
the O 3 column computed using po ps (Table 4; average col- 
umn O 3 difference from all sondes is +6.5 DU). 

Adding a typical O 3 climatology (e.g., McPeters and 
Labow, 2012) above-balloon burst allows calculation of to- 
tal O 3 column abundance for both the original and pressure- 
corrected ozonesonde profiles. In this case, offsets are re- 
duced to within a few DU (Fig. 9a). Note that sonde and/or 
satellite-based climatologies have become standard, replac- 
ing a constant mixing ratio assumption (McPeters et al., 
1997, 2007; Thompson et al., 2003; McPeters and Labow, 
2012; Morris et al., 2013). The constant mixing ratio as- 
sumption takes the O 3 mixing ratio at balloon burst and ex- 
trapolates that value to the top of the atmosphere to provide 
the above -burst residual column. Thus, O 3 MR errors such as 
those observed here will lead to significant errors in the resid- 
ual and total column O 3 if a constant mixing ratio method 
is used. The sonde-only O 3 column discrepancies brought 


about by the differences in the balloon burst altitudes be- 
tween the original and corrected pressure profiles are recon- 
ciled with the satellite climatological add-on above-balloon 
burst and comparison of the total column O 3 . The amount of 
total column offset is reduced to a mean offset within 3.8 DU 
for every radiosonde type with the above-burst addition (Ta- 
ble 4), signifying that both the O 3 MR error and altitude dif- 
ferences are contributing to total column discrepancies. 

A common practice within the ozonesonde community is 
to cut off total column O 3 integration at lOhPa (Thompson 
et al., 2003, 2007), rather than integrating the entire profile, 
and to apply a climatology such as that of McPeters and 
Labow (2012) to the remainder. This approach is employed 
and recommended for a variety of reasons including mitiga- 
tion of increasing pump efficiency uncertainties with altitude 
in the stratosphere (Johnson et al., 2002) and the reduced ac- 
curacy of the O 3 measurements above 10 hPa (Komhyr et al., 
1995b). The same technique was applied to the ozonesondes 
in this study to test if the sonde-only and total column O 3 off- 
set is reduced due to elimination of increasing O 3 MR errors 
routinely observed above lOhPa. 

The lOhPa cut off considerably reduces the differences 
between the uncorrected and pressure offset corrected sonde- 
only columns for most radiosonde types. Exceptions are the 
iMet-P launches, which in our data set rarely reached 10 hPa 
due to use of a smaller balloon (portions of the Houston 
and Idabel launches), and the RS92 profiles, which had lit- 
tle O 3 column error to begin with. Since the iMet-P launches 
rarely reached lOhPa, the entire balloon profile was inte- 
grated to the burst altitude, eliminating any effect this cut 
off would have had. With the exception of the iMet-P son- 
des, sonde-only column O 3 average differences are reduced 
from a maximum of 11.2 DU to within 3.5 DU (Table 4, 
Fig. 8 b) with a lOhPa cut off. Considering the total col- 
umn C> 3 with the lOhPa cut off and subsequent McPeters 
and Labow (2012) climatological add-on, the agreement be- 
tween the uncorrected and corrected pressure O 3 columns is 
further improved and most differences are essentially noise 
within the uncertainty of the total column integration from 
the ozonesonde. All radiosonde types agree to an average 
offset within —1.7 to +0.8 DU, with the poorest agreement 
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Fig. 6. Pressure offset at various altitudes versus eventual pressure offset at burst by radiosonde type. The various radiosonde types are 
identified by their respective colors. 
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dashed line marks the zero line for reference. 


from the iMet radiosonde 90th percentile of +3.5 DU (Ta- 
ble 4, Fig. 9b). 

Figure 10 shows analysis of an individual profile to un- 
derstand better the improved agreement in total column O 3 
after the pressure correction is implemented. It appears the 


standard lOhPa cut off may provide a serendipitous solu- 
tion to reconciling the differences between p and pops total 
and sonde-only column O 3 . The compensating effects of the 
pressure offset are viewed in terms of O 3 MR, P03 , and inte- 
grated sonde-only column with p and pops- Because lOhPa 
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Fig. 8. Ozonesonde-only column O 3 using the difference of 
columns calculated with p and pq pg, with integration to burst (a) 
and cut off at 10 hPa (b) compared to the 26 km pressure offset (p— 
POPS). The various radiosonde types are identified by their respec- 
tive colors. A few outliers were left from the figure for clarity. 

is above the po 3 maximum in the stratosphere, the discrep- 
ancies on either side of the O 3 peak routinely compensate 
for one another when sonde integration is truncated (i.e., the 
column differences below the O 3 peak are negative (positive) 
while above the peak they are positive (negative)). Integrat- 
ing to the burst altitude for those sondes that reach above 
lOhPa results in poorer agreement with altitude - the fur- 
ther above 10 hPa the sonde reaches before burst, the greater 
the column error becomes. Thus it appears that the lOhPa 
recommended limit for using the O 3 profile data results in a 
fortuitous minimization of the column errors caused by the 
pressure offsets and therefore our analysis argues in favor of 
the application of an O 3 climatology such as that used by 
McPeters and Labow (2012) above-balloon burst, with a cut 
off at 10 hPa if necessary. 

4 Summary and recommendations 

A total of 73 1 radiosondes were compared to quantify er- 
rors in radiosonde pressure sensor measurements relative to 


Fig. 9. As Fig. 8 except the McPeters and Labow (2012) above- 
burst O 3 climatology was added to each sonde from (a) burst or 
(b) 10 hPa/burst if greater than 10 hPa. 

pressure calculated from GPS measurements and to assess 
the impact these pressure offsets have on O 3 MR and column 
O 3 measurements. The pressure offset was shown to detri- 
mentally affect O 3 measurements, particularly in the strato- 
sphere, where errors in O 3 MR frequently exceed the labora- 
tory uncertainty of the ozonesonde measurements of ±5 % 
in the lower stratosphere (Komhyr et al., 1995b). The perfor- 
mance of Vaisala RS92 radiosondes was superior to RS80s 
and three series of iMet radiosondes, and was characterized 
by offsets of only ±0. 1-0.2 hPa at balloon burst, translating 
to O 3 MR errors generally within ±1-2 % at 26 km. The RS80 
and iMet-P radiosondes had the greatest 26 km average off- 
sets of — 1.01 and — 1.95hPa, respectively, translating to av- 
erage O 3 MR errors of —4.75 and —8.75 %. 

The differences between the radiosonde-measured and 
GPS-calculated pressures also introduced an altitude shift 
in the profile that must be considered for satellite valida- 
tion studies and column O 3 integration. The ozonesonde- 
only column exhibited a robust relationship with 26 km pres- 
sure offsets; sonde column differences between p and pavs- 
corrected profiles often exceeded ±10 DU, or ~ 3 % of the 
total column when offsets were beyond — l.OhPa at 26km. 
These column differences were reduced with the application 
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Fig. 10. Las Tablas, Panama RS80 sounding from 31 July 2007 
showing GPS (red) and radiosonde (blue) profiles of O 3 MR (a, left) 
and po 3 (b, right). The inset in (b) is integrated ozonesonde column 
showing compensating differences causing agreement in column O 3 
by 10 hPa. The 10 hPa cut off used prior to adding the McPeters and 
Labow (2012) O 3 climatology is marked by the black dashed line 
on all plots. The median 26 km pressure offset for this launch is 
-3.01 hPa. 

of the above-balloon burst O 3 climatology of McPeters and 
Labow (2012). When an integration cut olf of lOhPa was 
applied the agreement in total column O 3 between p and 
Pops profiles improved to within a few DU. The improved 
agreement between the uncorrected and corrected total O 3 
columns using a standard profile climatology and the lOhPa 
cut off argues for adopting this technique for column abun- 
dance estimates, especially with ozonesondes launched with- 
out GPS technology. Note that in the absence of GPS verifi- 
cation of the pressure profiles and O 3 MR, this cut off tech- 
nique only improves the resulting calculated column abun- 
dance and does not improve the accuracy of the O 3 profile 
shape or O 3 MR profile magnitude at the top of the profile. 

The ozonesonde community is currently in the process 
of homogenizing data (Deshler, 2012), seeking the highest 
accuracy trends and measurements, particularly at altitudes 
where satellite validation plays a vital role, from a global 
data set spanning dozens of stations and up to 40 yr of mea- 
surements. The homogenization process will take into ac- 
count sources of discrepancies and biases between different 
ozonesonde manufacturers, potassium iodide sensing solu- 
tion strengths, and pump efficiency corrections. The pres- 
sure offset introduces an additional source of error (often sig- 
nificant) that is independent of the ozonesonde partial pres- 
sure measurement, and an error that is not constant from one 
flight to the next, either with altitude or within a specific ra- 
diosonde manufacturer/type. It is anticipated that the analy- 
ses here will contribute to pressure corrections required as 
part of the ozonesonde data reprocessing. 



Fig. Al. Altitude differences between pressure and GPS altitude 
with GPS altitude based on magnitude of pressure offset. Pressure 
offsets from —3 hPa to 3 hPa in increments of 0.25 hPa were plotted. 
This calculation assumes a scale height of ~ 7000 m (T v = 240 K, 
g = gp = 9.80665ms -2 ), and was calculated from POPS — 1013 
to 5 hPa. 

The results of this study suggest the following recommen- 
dations regarding the pressure offset: 

1. Whenever possible, ozonesondes should always be 
launched with a GPS-enabled radiosonde to ensure an 
accurate O 3 MR magnitude and profile shape. 

2. Pressure-dependent variables should be recalculated 
using pgps, especially when pressure offsets exceed 
±1.0 hPa or ±5 % of the total atmospheric pressure at 
26 km/20 hPa. 

3. An above-burst climatology such as that used by 
McPeters and Labow (2012) should be applied us- 
ing a lOhPa cut off (if applicable), particularly with 
ozonesondes launched prior to the GPS era for column 
abundance observations. 

Appendix A 

The effect the pressure offset has on the difference between 
radiosonde-reported geopotential altitude and GPS altitude 
is presented in Fig. Al. Using standard gravity. g 0 , it is 
seen how a pressure offset of ±1.0 hPa (frequently observed 
in this study) can lead to an altitude discrepancy of ± 1 . 0 - 
1.5 km, having implications for column O 3 and a shifting of 
the O 3 profile shape. 

The pressure offset ( p-povs) and O3MR offset ([C>3mr(/?)- 
C>3mr(gps)]/03mr(gps)) by the launch site are shown in 
Figs. A2 and A3. As mentioned in Sect. 3.2, similar offset 
groupings are observed in the campaign-based launches from 
Porterville, CA, Las Tablas, Panama and at Idabel and Flous- 
ton which launched iMet-P sondes in the course of 10 months 
in 2012-2013. 
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Fig. A2. Pressure offset ( p-PGPS ) by launch site. A red dashed line marks the zero line for reference. 
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Fig. A3. Percent O 3 MR offset ([O 3 P-O 3 GP 5 ;]/ O 3 GP 51 ) by launch site. A red dashed line marks the zero line for reference. 
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